Dark Energy Constraints from the Cosmic Age and Supernova 



Bo Feng 1 , Xiulian Wang 2 , and Xinmin Zhang 1 
1 Institute of High Energy Physics, Chinese Academy of Science, 
P.O. Box 918-4, Beijing 100039, P. R. China and 



O 
O 

(N 
O 

Q 

> 

(N 
O 

o 

or 

6 
U 

in 
c3 



Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100080, P. R. China. 

Using the low limit of cosmic ages from globular cluster and the white dwarfs: to > 12Gyr, 
together with recent new high redshift supernova observations from the HST/GOODS program and 
previous supernova data, we give a considerable estimation of the equation of state for dark energy, 
with uniform priors as weak as 0.2 < fi m < 0.4 or 0.1 < n m h 2 < 0.16. We find cosmic age limit plays 
a significant role in lowering the upper bound on the variation amplitude of dark energy equation 
of state. We propose in this paper a new scenario of dark energy dubbed Quintom, which gives 
rise to the equation of state larger than —1 in the past and less than —1 today, satisfying current 
observations. In addition we've also considered the implications of recent X-ray gas mass fraction 
data on dark energy, which favors a negative running of the equation of state. 



Age limits of our universe are among the earliest mo- 
tivations for the existence of the mysterious dark energy. 
Namely, observations of the earliest galaxies could set 
a low limit on the age of the universe. In 1998, two 
groups 0, independently showed the accelerating ex- 
pansion of our universe basing on Type la Supernova 
(SNe la) observations of the redshift-distance relations. 
The recently released first year WMAP data 0] support 
strongly the concordance model with dark energy taking 
part of ~ 2/3. The most recent discovery of 16 SNe la 
4] with the Hubble Space Telescope during the GOODS 
ACS Treasury survey, together with former SNe la data 
alone could provide a strong hint for the existence of dark 
energy. Riess et al. 3 provided evidence at > 99% for the 
existence of a transition from deceleration to acceleration 
using supernova data alone. 

Despite our current theoretical ambiguity for the na- 
ture of dark energy, the prosperous observational data 
(e.g. supernova, CMB and large scale structure data and 
so on ) have opened a robust window for testing the re- 
cent and even early behavior of dark energy using some 
simple parameterization for its equation of state (e.g., 
Ref. Jp ) or even reconstruction of its recent density 
H S W Both recent WMAP fit and more recent fit by 
Riess et al. find the behavior of dark energy is to great 
extent in consistency with a cosmological constant. In 
particular when the equation of state is not restricted 
to be a constant, the fit to observational data improves 
dramatically 0, Hfl HE 113 • Huterer and Cooray [l(| pro- 
duced uncorrelated and nearly model-independent band 
power estimates (basing on the principal component 
analysis^!) of the equation of state of dark energy and 
its density as a function of redshift, by fitting to the re- 
cent SNe la data they found marginal (2-cr) evidence for 
W(z) < — 1 at z < 0.2, which is consistent with other 
results in the literature^, @, GjJ EJ El EE G3 

The recent fit to first year WMAP and other CMB 
data, SDSS and 172 SNe la data [H| by Tegmark et al 
E3 provided the most complete and up-to-date fit. Al- 
though SNe la data accumulated more after that, Ref. 
[T^| should still be a very profitable benchmark for cur- 
rent fit of the observables. However, when considering 



the behavior of dark energy alone, one has to do more 
since Ref. E3 only dealt with constant equation of state 
before the recent release of 16 more SNe la data by 
Ref. 4j. In fact a complete fit to full observational data 
still remains impossible provided one wants to recon- 
struct the full behavior of dark energy, despite the using 
of the most efficient Markov Chain Monte Carlo(MCMC) 
method. Under such circumstance Wang et al. and Riess 
et al. fitted dark energy to SNe la data, 2dfpfj| linear 
growth factor and a parameter related (up to a constant) 
to the angular size distance to the last scattering surface. 
In fact even the angular size distance is model dependent, 
as can be seen from Ref. 3] it differs for the six-parameter 
vanilla model and when an additional parameter a (run- 
ning of the spectral index) is added. Regarding the con- 
straint from the cosmic age, Krauss[2lJ used the WMAP 
fitted value for seven parameters: to — 13.7±0.2Gyr, to- 
gether with la HST [22| bound and assuming some spe- 
cific relation between f2 m and h, he got a lower bound on 
constant equation of state: W > —1.22, which is in strik- 
ingly agreement with WMAP result. Generally speaking, 
age limit can give an upper limit rather than lower limit, 
as shown by Cepa[23| . 

The low limit to the cosmic age can be directly ob- 
tained from dating the oldest stellar populations. Glob- 
ular clusters (GC) in the Milky Way are excellent lab- 
oratory for constraining cosmic ages. Carretta et al. 
|24| gave the best estimate for the age of GCs to be 
Age=12.9 ± 2.9Gyr at 95% level. The limit for age of 
GCs is around 11-16 GyrQ. White dwarf dating pro- 
vides a good approach to the main sequence turn-off. 
Richer et aZ.j25| and Hansen et a/.[2(| found an age of 
12.7 ± 0.7Gyr at 2a level using the white dwarf cooling 
sequence method. For a full review of cosmic age limit see 
Ref. 3]. The low limit to cosmic age serves as the 'anti- 
smoking gun' in excluding models which lead to shorter 
age. In this paper we use to > 12.0Gyr as the bound on 
cosmic age. Other constraints we use are only uniformly 
in range 0.2 < Q m < 0.4 or 0.1 < VL m h 2 < 0.16. As 
can be seen from Ref. 0] > such constraint is much looser 
than the six parameters + W set and comparable to the 
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nine-parameter- vary set 1 . We assume this constraint 
to be reasonable with one additional parameter added 
below: the variation of equation of state. We find for the 
linear parameterization of W age constraint can shrink 
the upper bound on W from ~15 to ^5 when using SNe 
la alone and from ^5 to ^2 when considering above pri- 
ors on f2 m or il m h 2 . While for the model introduced by 
LindeQ, the upper bound on W a can shrink from ^TO 
to ^5 when with priors on Q m or fl m h 2 . 
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FIG. 1: Age and SNe la constraints on ACDM cosmology. 
Left panel: The red area is excluded by to > 12Gyr and 
the blue area is excluded by the assumption for to < 20Gyr 
for ACDM model. The area between the two solid lines is 
allowed by la HST limit. Right panel: 2a SNe la limit on 
ACDM model. The dashed line corresponds to the la limit 
and the dot inside denotes the best fit value. The navy area 
is allowed by age constraint 12Gyr < to < 20 Gyr. 



The cosmic age can be written as 

t -H-i r dz 

° ° Jo (1 + *)£(*) ' 

where Q 

e(z) = [n m (i + z) 3 + (i - n m )x( z )] 1/2 

and X(z) = p x (z)/p x (0). 



(1) 



(2) 



1 Our prior on Q m is also consistent with the median statistics 
study on mass density by Chen and Ratral23 . fo r more investi- 
gations on the effects of priors see Refs.Hl l28l.l29l . 



Firstly we delineate the effect of age limit in ACDM 
cosmology. In the full paper we assume a flat space, i.e. 
fit = 0. In the left panel of Fig. 1 we vary Vl m from 
to 1 and the Hubble parameter h from to 1.4. The red 
area is excluded by to > 12.0Gyr, the area between the 
two black solid lines is given by the la HST limit. If we 
conservatively assume in ACDM cosmology the cosmic 
age is no more than 20 Gyr, the area with blue color will 
be excluded. It can give a rough estimate for current 
fraction of matter in the universe, although much looser 
than SNe la constraint as shown in the right panel. The 
supernova data we use is the "gold" set of 157 SNe la 
published by Riess et al. in Q . In the below we constrain 
the Hubble parameter to be uniformly in 3a HST region: 
0.51 <h< 0.93. 

In the detailed discussions below we consider two type 
of parameterizations for the equation of the state of the 
dark energy, 



W(z) = W + W'z 7 



(taken as Model A) so that 



— n -3(l + W-W ) 3W' (a,- 1 -1) 



X = a 



(3) 



(4) 



where a is the cosmic scale factor. The other form was 
proposed by Refs. |5j, |3CJ (taken as Model B): 



W[z) = Wi + W a z/{1 + z) 



which leads to 



X 



a -3(l+W 1 +W a ) e 3W a (a-l) 



(5) 



(G) 



Both models(as well as other models such as firstly pro- 
posed in [LUp make good approximations to probe the 
behavior of dark energy around the present epoch, while 
the former model leads to poor parameterization at very 
large redshift. But as argued by Riess et al. Q this is 
acceptable for showing the late behavior of dark energy. 

SNe la data alone proves to be a weak constraint 
on above models.as shown in the right panel of Fig. 2 
and also in Ref . Ml where the authors used flux aver- 
aging method [Til l3l|. In the left panel of Fig. 2, the 
region on up right corner is fully excluded by the age 
limit(<o > 12Gyr). The role of age limit can be easily seen 
from Eqs.(l-6). In Model A larger W' leads to larger X 
in early epochs, hence corresponds to smaller ages which 
can be directly constrained. Similar case works for Model 
B. Age constraint still works when adding the prior on 
fl m or fl m h 2 . In Figs. 3-4 we show the corresponding 
effects when adding different priors for SNe la on right 
panels and show the role of age correspondingly in the left 
panels. The dashed lines are the lcr regions and the small 
dots denote the best fit parameters. Up right regions of 
the left panels are excluded by the cosmic age limit. We 
can see the age limit reduces significantly the upper re- 
gions and consequently changes the best fit values of the 
model parameters. One can also find from Fig. 2 and 
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FIG. 2: Right panel: 2<r SNe la limit alone on Model A dark 
energy. Left panel: 2a SNe la limit and age limit (to > 12Gyr) 
on Model A dark energy. The dots inside the two panels show the 
best fit parameters. 

Fig. 3 the role of priors in deriving w, although the ef- 
fect of age limit does not change, the priors on fi ro has 
shrinked the allowed parameter space on w. Generically 
the results on w also depends on different parameteriza- 
tions, as shown in our paper and also in the literature(e.<7. 
Ref. However this would not change the picture 

that cosmic age can put an additional constraint on the 
equation of state w. 

As mentioned above, the present data seem to favor 
an evolving dark energy with the equation of state being 
below —1 around present epoch evolved from W > — 1 
in the past. This can also be seen from the best fit pa- 
rameters of our Figs. 2-5. If this result is confirmed in 
the future, it has important implications for the theory of 
dark energy. Firstly, the cosmological constant as a can- 
didate for dark energy will be excluded and dark energy 
must be dynamical. Secondly, the simple dynamical dark 
energy models considered vastly in the literature like the 

quintessence j^E or the phantom^, HE HE can 
not be satisfied either. 

In the quintessence model, the energy density and the 
pressure for the quintessence field are 

P=\Q 2 + V{Q) , p=\Q 2 -V{Q) . (7) 

So, its equation of state W — p/p is in the range 
-1 < W < 1 for V(Q) > 0. However, for the phan- 
tom ~|34| which has the opposite sign of the kinetic term 
compared with the quintessence in the Lagrangian (we 
use the convention (+, — , — , — ) for the sign of the met- 
ric), 

C = ~d fi Qd»Q-V(Q) , (8) 

the equation of state W = ( - \ Q 2 - V) / ( - \ Q 2 + V) is lo- 
cated in the range of W < — 1. Neither the quintessence 
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FIG. 3: Age and SNe limits on Model A with different priors as 
noted inside. The dots inside the la dashed lines denote the best 
fit parameters. 

nor the phantom alone can fulfill the transition from 
W > — 1 to W < — 1 and vice versa. 2 But at least a 
system containing two fields, one being the quintessence 
with the other being the phantom field, can do this job. 
The combined effects will provide a scenario where at 
early time the quintessence dominates with W > — 1 and 
lately the phantom dominates with W less than —1, sat- 
isfying current observations. As an example, we consider 
a model: 

-F [exp(-— 0i) + exp(-— fa)] , (9) 
m p m p 

where fa and fa stand for the quintessence and phantom. 



2 Although the k-essenceQ like models can have W < — 1|23, it 
has been proved later by Ref . |4flj to be difficult to get W across 
-1 during evolution. 
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FIG. 4: The same as Fig.3 for Model B. 



In Fig. 0] we illustrate the evolution of the effective equa- 
tion of state of such a system with — ln(l + z). 

In general to realize the transition of W around — 1, 
one needs to consider models of dark energy with more 
complicated dynamics and interactions with gravity and 
matter. This class model of dark energy, which we dub 
"Quintom" , is different from the quintessence or phan- 
tom in the determination of the evolution and fate of 
the universe. Generically speaking, the phantom model 
has to be more fine tuned in the early epochs to serve 
as dark energy today, since its energy density increases 
with expansion of the universe. Meanwhile the Quintom 
model as illustrated in Fig. 5 can preserve the tracking 
behavior of quintessence |4ll |42| . where less fine tuning is 
needed. We will leave the detailed investi gati on of the 
Quintom models in a separated publication |43j ] , however 
will mention briefly two of the possibilities below in ad- 
dition to the one in Eq. (9). One will be the scalar field 
models with non- minimal coupling to the gravity 0, 
where the effective equation of the state can be arranged 
to change from above -1 to below -1 and vice versa. For 
a single scalar field coupled with gravity minimally, one 
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FIG. 5: The evolution of the effective equation of state of 
the double scalar fields given in Eq. The parameters are 
chosen as: V = 8.38 x 10" 126 mi A = 20. We set the initial 
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may consider a model with a non-canonical kinetic term 
with the following effective Lagrangian |43j : 



C 



lf(T)d^Qd»Q 



V(Q) , 



(10) 



where f(T) in the front of the kinetic term is a dimen- 
sionless function of the temperature or some other scalar 
fields. During the evolution of the universe when f(T) 
changes sign from positive to negative it gives rise to an 
realization of the interchanges between the quintessence 
and the phantom scenarios. 

Recently Allen et al. 46] have provided new obser- 
vational data basing on Chandra measurements of the 
X-ray gas mass fraction in 26 X-ray luminous galaxy 
clusters. Under the assumption that the X-ray gas mass 
fraction measured within ^2500 is constant with redshift 
the fgas data in the range 0.07 < z < 0.9 can be used 
directly to constrain cosmological models. We use their 
data and fit to Model A, we set the same gaussian prior 
on flbaryon as Ref. |46j meanwhile varying h uniformly in 
range 0.51 ~ 0.93 and < O m < 1. The x 2 value is 
defined as 

[/gas ( Z i) — /gas,i] \ 
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FIG. 6: 2— dimensional and 1— dimensional constraints on 
W' in light of the data f gas from Chandra observations from 
Ref. lirl with SNe la constraints also shown for comparison. 



ft h h 2 - 0.0214 V 



0.0020 



b -0.824V 



0.089 



(11) 



For a detailed description of the data and fitting see 

Ref. |4(j . We delineate the resulting 2— dimensional and 
1— dimensional plots in Fig. 6, together with previous 
SNe la results for comparison (here we do not include 
age constraint for simplicity). Very interestingly we get 
W < with the center value (W, W) = (-0.2, -4.875) 
when varying fine grids, which indicates the universe 
may not be accelerating today but in the very near 



past. However the smallest redshift of f gas data gives 
z = 0.077900^(1, this shows also some poor parame- 
terization of model A. As shown in our paper and also 
vastly in the literature the results of fitting depends some- 
what in the parameterizations |50j . However, in any case 
if nonzero W gets more favored with the accumulation of 
observational data, it gives strong implications for dark 
energy "metamorphosis" 14J. We also find although the 
two data sets give consistent results to model A, there 
seems to be some discrepancy. In Fig. 6, the contours 
do overlap in la regions, but only in a small area and 
the right panel is more distinctive: f gas favors a neg- 
ative W at more than 1.3a while a positive W' is fa- 
vored at around l.lcr with the prior on f2 m for SNe la. 
Basically the X-ray data probe the late behavior of the 
angular-diameter distance and SNe la probes the lumi- 
nosity distance. There seems to be some discrepancy 
between them and this may possibly be due to some new 

physicslHIHIH. 

In summary in this paper we consider the effect of cos- 
mic age and supernova limits on the variation of W . Our 
results show that age limit plays a significant role in low- 
ering the variation of amplitude on the equation of state. 
Current SNe la observation seems to favor a variation of 
W from > —1 in the recent past and < — 1 today. If 
such a result holds on with the accumulation of obser- 
vational data, this would be a great challenge to current 
cosmology. We give a simplest example of Quintom which 
can satisfy the current implications on the equation of 
state on dark energy, and discuss briefly the possibility 
of building Quintom models. 
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